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We present a detailed analysis of the micropipet experiments recently reported in J-B. Manneville 
, et al., Phys. Rev. Lett. 82, 4356-4359 (1999), including a derivation of the expected behaviour 

' of the membrane tension as a function of the areal strain in the case of an active membrane, i.e., 

^SJ I containing a nonequilibrium noise source. We give a general expression, which takes into account 

> I ■ the effect of active centers both directly on the membrane, and on the embedding fluid dynamics, 

' keeping track of the coupling between the density of active centers and the membrane curvature. 

The data of the micropipet experiments are well reproduced by the new expressions. In particular, 
we show that a natural choice of the parameters quantifying the strength of the active noise explains 
' both the large amplitude of the observed effects and its remarkable insensitivity to the active-center 

04 , density in the investigated range. 
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Biological membranes are made up of a complex mixture of lipids and proteins. The lipid molecules form a bilayer 
structure which separates the cytoplasm of the cell from the outside. In addition to this structural role, the membrane 
also participates in a number of the living cell functions mostly performed by proteins embedded inside the lipid 
^ ' bilayer, such as solute transport via ion channels or pumps, cell locomotion and adhesion, membrane transport 
through exocytic and endocytic pathways, signal transduction... Consequently, from the statistical physics point 
of view, biological membranes are strongly out of thermodynamic equilibrium, whereas most studies on membranes 
reported in the physics literature have been done at thermodynamic equilibrium In order to achieve a more 
^vq complete physical description of biological membranes, this nonequilibrium aspect clearly has to be included. The 
^vq field of membrane shape fluctuations is a good test case in which to examine the relevance of nonequilibrium effects. 
At thermodynamic equilibrium, the membrane shape fluctuates because of thermal noise, i.e., the Brownian motion 
of the bilayer. Such membrane will be called a 'passive' membrane in this paper. If a nonequilibrium noise source 
[ is superimposed to thermal noise, due for instance to the activity of membrane proteins, then the membrane is no 
\^ • longer at thermodynamic equilibrium. In that case, the membrane will be called 'active'. 

Recently, micropipet experiments on fluctuating giant vesicles containing bacteriorhodopsin (BR) reconstituted 
in the lipid bilayer have shown that the light-driven proton pumping activity of BR induces an amplification of 
the membrane shape fluctuations In these experiments, a larger excess area could be pulled out by micropipet 
^ aspiration when the proteins were activated. The results were qualitatively interpreted in terms of an increase 
'"^ of the effective membrane temperature and were not directly compared to theoretical predictions. In the present 
^ [ article, we give details about the experimental procedure (section ||) and develop a theoretical framework to analyze 
O ■ quantitatively the micropipet experiments (sections III and |^) . According to theory, a qualitatively new fluctuation 



spectrum is expected in the presence of a nonequilibrium noise source |^^. These earlier theories introduced the 
^ : nonequilibrium activity in the membrane conformation equation only. This restrictive choice was made because 
' the nonequilibrium force-density arising from the active proteins, when included in the Stokes equation for the 
, solvent velocity field, altered the membrane fluctuation spectrum through terms which were subdominant at small 
• wavenumber. 

We show here that those nominally subdominant terms provide the most important contribution in the experimen- 
tally relevant range. This unexpected behavior is due to the very small value of the permeation coefficient. With this 



implementation, theory and experiment are brought in agreement (section IV). Even the absence of sensivity of the 



experimental data on the active center density appears as a natural consequence of the developed theory. 
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II. EXPERIMENTAL PROCEDURE AND RESULTS 



A. Bacteriorhodopsin 



The bacteriorhodopsin (BR) is a 27 kDa protein purified from the so-cahed purple membrane of the halophihc 
bacteria Halobacterium salinarum Its structure is known at the atomic level with high resolution |^. The BR 
absorption spectrum shows a maximum in the green- yellow wavelength around 566 nm (Fig. |l|) . 




FIG. 1. Absorption spectrum of the purple membrane suspensions used in the experiments (solid line) superimposed with 
the absorption spectra of the high-pass green-yellow filter (dashed line) and the high-pass red filter (dotted line). 

When illuminated with green-yellow light, proton pumping is activated through a photocycle involving several 
distinct photointermediates jl^. The total duration of the photocycle is r ~ 5 ms. Structural changes of the BR 
during the photocycle have been investigated to elucidate the translocation pathway of the proton across the protein. 
The pumping mechanism has been recently completely elucidated, so that BR is to date the best understood ion pump 
^d| . The proton pumping activity has been extensively studied in reconstituted systems, mostly in large unilamellar 
vesicles (0.1 — 1 /im in diameter) JT^ . 



B. Giant vesicle formation 

In referenc e fl , the electroformation technique of giant unilamellar vesicles (10 — 100 /im in diameter) jl^, modified 
according to |14[| for BR incorporation in the lipid bilayer, was used to grow giant vesicles from a mixed lipid/protcin 
dried film. The phospholipid EPC (Egg Phosphatidylcholine; Avanti Polar Lipids, Alabaster, Alabama, USA) is 
a mixture of lipids with different chain lengths and degrees of saturation and is known to be adequate for BR 
incorporation [|5|. EPC (0.5 mg/ml) was first resuspended in diethyl ether. Concentrated BR (18 mg/ml) was then 
added at a molar ratio of 80 lipid molecules per BR molecule. The mixture was sonicated at °C for 30 seconds 
and a few microliters were deposited on ITO (Indium Tin Oxide) treated glass slides at 4°C. The protein/lipid film 
was dried under vacuum overnight. A vesicle electroformation chamber was formed by assembling and sealing with 
wax (Sigillum wax; Vitrex, Copenhagen, Denmark) two ITO slides separated by 1 mm Teflon spacers. The fllm was 
hydrated by injecting a 50 mM sucrose solution in the chamber. An electric field (1.5 V AC) was applied across the 
chamber by connecting the ITO slides to copper electrodes. Giant vesicles were obtained in about two hours and 
transferred in a micromanipulation chamber filled with 50 mM glucose to enhance the optical contrast between the 
inside and the outside of the vesicles. Sodium azide (1 mM) was first added to the sucrose and glucose solutions 
to avoid bacteria proliferation. In some experiments, respectively 16 and 25 % (w/w) glycerol was added to both 
the internal and external solutions in order to increase the viscosity to respectively 1.577^, and 21]^, where rj^ is the 
viscosity of water. 

BR incorporation was checked by fluorescent labelling of BR with FITC (Fluorescein Isothiocyanate, F-143; Molec- 
ular Probes, Eugene, Oregon, USA) following a published protocol |l^. Excitation of the FITC was performed at 
488 nm with an argon laser (Spectra Physics, Les Ulis, France) through the epi-illumination pathway of an inverted 
microscope (Axiovert 135; Zeiss, Oberkochen, Germany). The fluorescence images of the vesicles were acquired by a 
low light level SIT (Silicon Intensified Target) camera (LII4036; Lhesa, France) (see Fig. ^). The fiuorescence inten- 
sity If — (Ives — hgdl/hgdi where I^es is the fluorescence intensity of the vesicle contour and hgd is the background 
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intensity, was measured by computer image analysis using a CH — h custom software running on a Pentium 200-based 
computer with a Meteor frame grabber (Matrox, Rungis, France). 
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FIG. 2. Experimental set-up designed for the micropipet experiments. 



The BR was activated by illumination through a high-pass (500 nm) green-yellow filter located in the trans- 
illumination pathway (see Figs. |^ and H), to avoid the non-pumping branched photocycle initiated if the M interme- 
diate absorbs at 440 nm |T^. To image the vesicles without activating the BR, a high-pass (650 nm) red filter was 
substituted for the green-yellow filter (see Figs. and ||). The illumination power was in the same range as that 
known to fully activate BR reconstituted in large unilamellar vesicles (lO'^W/m^) jf^]. To correct for the different 
bandwidths of the green- yellow and red filters, the trans-illumination light focused on the specimen plane was adjusted 
to about lOOmW/cm^ in all the experiments. The sample was illuminated for at least 15 minutes before starting an 
experiment, so that the BR was always in its light-adapted form fT^ . 

It has been shown that the reverse phase evaporation technique used to incorporate BR in large unilamellar vesicles 
(typically 200 nm in diameter) results in an asymmetrical orientation of the BR molecules across the lipid bilayer 
. Consequently, for these vesicles, a proton gradient builds up across the lipid membrane upon activation, which 
inhibits BR pumping activity. However, since the electroformation technique is symmetrical, we do not expect any 
asymmetry in the BR orientation, and thus we do not expect inhibition of the pumping activity. To be on the 
safe side, we have performed additional experiments which were designed to cancel any proton gradient according to 
the following procedure. A classical way of suppressing the inhibitory proton gradient, without incorporating any 
additional active molecule in the membrane, is to add KCl (potassium chloride) to the solution, since protons can 
then codiffuse passively through the bilayer in the form of HCl, and since chloride ions can diffuse inside the vesicle 
to ensure electroneutrality. In reference KCl was added to both the internal and external solutions up to 2 mM, 
a concentration above which electroformation of giant vesicles fails, in order to get rid of a possible proton gradient. 
Our results proved to be insensitive to the addition of KCl . 
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C. Micropipet experiments 



The micropipet technique developed by Evans and coworkers [[l8[ allows a quantification of the excess surface area 
stored in the membrane fluctuations by pulling it out with a micropipet aspiration: when a pressure difference is 
applied, the membrane is put under tension and sucked inside the pipet. The experimental set-up was built on 
an inverted microscope equipped with a 40x objective (N.A. 0.75 air Ph2 Plan Neofluoar, Zeiss) for epi-fluorescence, 
phase contrast and differential interference contrast (DIG) microscopy. The transmission phase contrast or DIG images 
were recorded by a GCD (Charge Coupled Device) camera (Sony, Paris, France). The sample cell was temperature 
controlled at 15 °C by a water flow to limit evaporation of the solution (see Fig. ^, bottom). 




FIG. 3. Typical micropipet experiment. The pressure difference AP is the difference between the pressure outside and 
inside the pipet. The intrusion length L is the length of membrane aspirated inside the pipet when a pressure difference AP 
is applied. The bar represents 10 /im. 

Glass micropipets were pulled from 1 mm outer diameter borosilicate capillaries (GC lOOT-10; Phymep, Paris, 
France) with a micropipet-puUer (P-97; Sutter Instruments Co., Novato, California, USA). The micropipet tip was 
cut on a microforge to obtain diameters up to 5 — 10 /im. The pipets were treated with BSA (bovine serum albumin, 
1 %) for 30 minutes to prevent adhesion of the lipid membrane to the glass pipet walls. A pipet holder was mounted 
on a three-dimensional piezo micromanipulator stage (Physik Instrumente, Waldbron, Germany) in order to control 
the position of the pipet tip within 0.1 /im accuracy. The pressure difference AP between the outside and the inside 
of the pipet was measured by a liquid-liquid pressure transducer (DP103-20; Validyne, SEI3D, France) with 0.01 
Pa accuracy. The pressure is imposed by a water height difference between two water filled tanks equipped with 
micrometric displacements (see Fig. ||, bottom). The absence of any air bubble in the water circuit running from the 
tanks to the micropipet is crucial and was checked before each experiment. The relationship between the pressure 
difference AF and the imposed membrane tension a directly derives from Laplace's law ]l^ : 



2(1 Rpip I Rves) 



where Rpip is the pipet radius and Rves is the vesicle radius, both measured directly from the DIG image (see Fig. |^). 

The excess area stored in the membrane shape fiuctuations a is defined as a — [A — Ap)/A, where A is the actual 
area of the fluctuating membrane and Ap is the area projected on the mean plane of the membrane. During a 
micropipet experiment, the excess area decreases as the membrane undulations are pulled out by the increasing 
pressure difference. For a given AP, an intrusion length L is aspirated inside the pipet (see Fig. ^ . A reference state 
(APq, Lq) is defined as the lowest suction pressure that can be applied in the experiment to aspirate the fluctuating 
vesicle inside the micropipet ||l^. The variation of the excess area as compared to this reference state, the so-called 
areal strain Aa = — a, follows from geometrical considerations. To first order in AL = L — Lq, we have |l|]: 

A (Rpip/ Rves) (Rpip/ Rves) a 7- 

Aa = ao - a = — AL 

^J^pip 

The increase of the intrusion length AL was measured by image analysis with pixel accuracy, i.e. 0.2 /im with the 
40x objective. 

The excess area can be expressed using the local displacement u of the membrane around its mean plane p9| : 
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In the low q regime the convergence of the integral is guaranteed by a crossover from a curvature dominated regime 
(u(q^)^) — kT / Kq\ to a tension dominated regime {u{(\±)'^) = kT/aq\ for a passive membrane, where k is the 
Boltzmann constant, T is the absolute temperature and k, is the bending modulus of the membrane. The upper limit 
is qmax = 2.1: /a, where a is a microscopic length. In the entropic regime, i.e. at low tension, inserting the fluctuation 
spectrum of an equilibrium membrane gives the dependence of the excess area a as a function of the membrane tension 
a p8| , p^ : a — (kT/Snfi) \n{cst/a), where est is an integration constant. The areal strain Aa = ao — a is thus: 

kT a 

Aa — aQ — a — - — In — . (2) 

SVTK (To 

For a passive membrane, the linear relationship between the logarithm of the tension In a and the areal strain Aa 
allows the determination of the bending modulus k. We will give in section |^ a similar relation relevant to the active 
case. 



D. Essential results 

The results reported in and duplicated in Figure |4| show that when the vesicles are illuminated with green- yellow 
light, the slope of the logarithm of tension versus the areal strain is smaller than in the case where the vesicles are 
illuminated with red light. This indicates that the excess area is larger when the BR is illuminated with green- yellow 
light, and consequently that BR activity induces an amplification of the membrane shape fluctuations. The quality 
of the fit suggests that one can describe the effect of the BR activity in terms of an effective temperature T'^^^ ~ 2T. 
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FIG. 4. Variation of the logarithm of the tension a versus the areal strain Aa for the same vesicle containing BR, 
alternatively passive and active. 

An other important feature of the experiment concerns the dependence of T^^ ^ on BR concentration. Fig. ^ 
shows that in a concentration range that we estimate between approximatively 10^^ and lO^^BR/m^ the effective 
temperature is essentially independent on BR concentration. This may look surprising since the BR activity is the 
driving force. 
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FIG. 5. Variation of the effective temperature T'^^-^ of the active membrane as a function of the fluorescence intensity If, 
thus of the BR concentration. 



Before developing the interpretation of these results, we must first guarantee their reliability, i.e. that it is an effect 
related to the out-of-equilibrium pumping activity of BR and nothing else. The control experiments with pure lipid 
vesicles obviously exclude a role of the lipids themselves. For these, we find for both green-yellow and red illumination 
the expected kT/K ~ 0.1 value [Q. Using simple estimates, we can also rule out the possibility of any thermally 
induced artefact due to the larger absorption of light by BR in the green-yellow wavelength. Assuming that one BR 
molecule absorbs one photon each r ~ 5 ms, the total stationary flux (total energy received per unit area of membrane 
and per unit of time) is: 

W = — X p 

T 

where h is the Planck constant, v is the photon frequency and p is the mean BR density. In a pessimistic estimate, we 
assume that this total flux W is dissipated via conduction in the surrounding water. The sample cell is temperature 
controlled by a cold water circulation and we assume that a temperature gradient arises from the BR heating between 
the membrane and the sample cell wall. This gradient extends over a typical length L = 1 mm which is the size of 
the sample cell. With C — 4.18 x 10^ J.m^^. °K~^ the heat capacity of water and Kt = 1-5 x 10^''m^/s the heat 
diffusion coefficient, we have: 

CKr^^W^^xp 
L T 



and 



hv 

AT = — X p X 



T CKt 

This yields a temperature increase of AT — 2 x 10^'^°K, five orders of magnitude smaller than the reported increase 
in effective temperature, which can not account for the observed effect Direct heating is just totally inefficient 
(note also that direct heating of water is clearly ruled out by the experiments on pure phospholipidic vesicles). 

Most importantly, the experiments with glycerol prove that the observed effect is of nonequilibrium origin. The 
addition of glycerol modifies dynamic parameters such as the solvent viscosity rj, the permeation coefficient Ap and 
the active force Fa. At thermodynamic equilibrium, such parameters cannot play a role in the fluctuation spectrum 
as imposed by the fluctuation-dissipation theorem. For an active membrane however, these parameters play a role as 
can be seen from references |Q] or 1^. The addition of glycerol increases the solvent viscosity while it decreases its 
permeation coefficient. The results given in reference report a lower increase in the effective temperature when 
16% and 25% (w/w) glycerol is added, clearly reveahng the out of equihbrium nature of the effect. This result is 
qualitatively consistent with the observation that the BR pumping activity is diminished upon addition of glycerol 
due to an increase in the lifetime of the M intermediate Finally, the fact that the observed effect does not 

depend on the measuring sequence (red light experiment or green light experiment first) rules out a potential role 
of the conformational change between the light-adapted and the dark-adapted states This is also consistent 

with the result that the ratio kT/ k in the passive case is the same as that of the pure phospholipidic vesicles. The 
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renormalisation of the bending rigidity by the BR is not measurable. All these observations give strong support to 
the assertion that the effect is indeed due to the proton pumping activity. 

The use of an effective temperature to qualitatively interpret the results according to equation is justified by 
the good quality of the linear fits of the micropipet experiments performed in j^. However, we need to develop a 
complete theoretical analysis to understand all these experimental features quantitatively. 

III. THEORY 

As in references we consider situations in which a membrane under tension is subjected to random forces 

of two different origins. These arise (i) from thermal agitation, i.e., the Brownian motion of the membrane, simply 
because membrane and solvent have a thermodynamic temperature, and (ii) from 'biological' activity such as proton 
pumping of the BR. The membrane equation of motion can be written to lowest order: 

\p\-^{r^,t) - K(r±,t)) = 6Pir^,t) - SU{r^,t) + Fai>ir^,t)+F;^pAMr±,t) + fpir^,t) (3) 

In this expression, u(rj_, t) is the membrane displacement at point rj_ = (x, y) with respect to a (x, y) reference plane 
orthogonal to z, the average membrane normal and Aj_ is the Laplacian in the xy plane. Vz{r±,t) is the fluid velocity, 
in the normal direction at the membrane surface; Ap is the membrane permeation coefficient. 6P{r±,t) — P{r±,z ~ 
0"*", t) — P(rj^, z = 0~, is the pressure difference across the membrane and ^n(rj^, t) ~ n(rj^, z = 0+, t) — n(r^, z = 
0~, t) the osmotic pressure difference. This osmotic pressure difference results from the proton pumping activity: for 
each BR cycle one proton is transferred across the membrane. In principle the calculation of SIl cannot be achieved 
without solving all dynamical equations of the problem. However, considering the convective term of theproton flux 
as a second order correction allows to evaluate SIl separately. We postpone this derivation to appendix ^ fp is the 
brownian noise acting on the membrane corresponding to the dissipation of energy in the permeation process and 
satisfies: 

(/p(ri,t)) =0 

(/p(r±,<)/p(rl,t')) - 2kT\p'5{r^-r'^)5{t-t') 

The term Faip{r±, t) results from the BR activity. More precisely ij}{Y±,t) — p^(rx, t) — p^{r±,t) is the local difference 
between the density {r± , t) of BR molecules transferring protons in the direction z (up) and the density p^ {r± , t) 
of BR molecules transfering protons in the -z direction (down) . Fa is the average elementary force transmitted to the 
membrane by a steady proton transfer. The flip-flop of BR is expected to be much slower than that of phospholipids 
and thus p^ (r±,t) and p^(r_L, t) can be considered as separately conserved quantities. As explained in section |l|, with 
our experimental conditions, the probability of inserting a BR molecule into the phospholipid membrane does not 
depend on the pumping direction, thus {p^) = (p^) and = 0. XpFa can be understood as measuring the average 
volume transferred through the membrane per BR and per unit time. Thus, XpFa^p(r±,t), which we will refer to as 
the 'active-permeative' term, can also be understood as the local volume transferred through the membrane per unit 
area and per unit time due to the pumping imbalance between up and down BR molecules. Terms corresponding 
to the stochastic nature of the pumping activity have been omitted since they have been shown to lead to smaller 
effects than those due to the collective iIj fluctuations |^,||. As discussed in Q, the fourth term describes the fact 
that pumping may work better when the membrane is curved with a given sign: F^ measures this sensitivity per 
pump and p — p'^{r±,t) + p^{r±,t). There is experimental evidence for such a coupling in the functioning of certain 
ion-channels called TRAAK and TREK H. 

We further need equations for the fluid and for the BR density dynamics. Navier-Stokes equations have to be 
implemented in two ways: one first has to keep track of the Laplace force exerted by the membrane on the fluid, and 
this can be done in the usual way 

Secondly, each BR has a small but finite spatial extent. Its activity will disturb the ambient solvent in the form of 
a distribution of force densities in its vicinity. Since no external force source is present, the total force must vanish, 
but its first moment will in general be present. For convenience, we adopt the simplest set of force-densities consistent 
with this requirement: a pair of oppositely directed point forces, separated by a distance of order the size of a BR 
molecule. This implies a dipolar force density Fa [S{z — w^) — 6{z + w^)] tp{r±, t) in the Stokes equations (this term 
will be called the 'active-hydrodynamic' term), and are lengths of the order of the BR size; their values are a 
priori unequal since the BR, or any molecule with unidirectional activity, is not up-down symmetric. Similarly, the 
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term [S{z — w'^) — 6{z + w^] pA±u{r±,t) describing the sensitivity of the force dipole to curvature should be kept. 
Thus, in the low Reynolds number regime appropriate to these experiments, we can write: 

= -VP(r,t)-|f(r^,t)5(z)z 

+Fa [S(z - w^) - S{z + w^)] tjj{r, t)z 

+i^^A_LM [S{z -w^)-S{z + w^)] p(r^, <)z 

+77AV(r,t)+f;,(r,i), 

where r refers to the three-dimensional position vector, and has the same meaning as in (|). P{r,t) is the 
three-dimensional pressure field, V(r, t) the three-dimensional fluid velocity field, F is the membrane free energy: 

F^IJ d^^^ [K{AMr±))^ + (T(Vu(ri))2 - 2E^{r^)A^u + x^'(r±)] 

K is the membrane 'bare' bending modulus, a the membrane tension, S a coefficient linking membrane curvature and 
BR imbalance, x the 'bare' susceptibility corresponding to that imbalance (for small enough densities x ~ ^T/p). 
Orders of magnitude will be given in the next section. The third term of equation (^), is the 'force dipole' density 
already described, the fifth and sixth are the usual viscous terms and associated forces: 

(f/.(r,t)) -0 

(A.(r, t)hj{v', t')) = 2kT7^ {-%V2 + d,dj}5{v - r')5{t - t') 

Last, we need a dynamical equation for the BR imbalance density ijj. Following we can write in the linear 
regime: 

with A = D/xi where D is the diffusion coefficient of the BR molecules in the membrane. This expression is valid for 
— 0, in the absence of fluctuation corrections, and the last term of (^ is a conserving noise, i.e., the divergence 
of a random current with 

(fv,(ri,t)) =0 

{U,{v^,t)U,{v^',t')) = 2AfcT% 6{v^ - v'^)5{t - t') 

In order to compare experiment and theory, we need to calculate the equal time correlation function 
(w(q_L, t)it(q'^, i)). We first eliminate Vz in (^) by solving for it from the Stokes equation (||) in Fourier space to 
get: 

— (qi,0 + T-^u(qi,t) = /3VXq±,t) + M (6) 

'di 

in which we have chosen the convention 



(qi,t) -i-T^ V^q^,*) = 7"(q±,*) + (7) 



/(q_L,i) = J .f{r_L,t) exp(iq^ ■r^)d^rj 



for Fourier transforms, and used 511 ~ 2kTaij}{(\^, t)/Dpq± where d is the proton pumping rate and Dp an effective 
proton diffusion coefficient (see appendix ^) . The parameters entering equations (||) and (0) are listed below: 

^u' = (AP + j^K<rql + nqi) + pXpF^ql - ^ql 



= Kxql = Dql F^ = F^ - 2kTa/Dpq^ 



P=XpFa-VaH^-^ql{Xp + J^) Va^Fa 

M = Ap/,(q^,t) + ^ / '-^dqz V'^ = Fi^ 

V = -«q_L fv)(q±,t) 
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Calculating the equal time {u{q±,t)u{q'^,t)) correlation function is straightforward although somewhat tedious; we 
find 



kT 

(u(q_L,t)w(-qi,t)) = 



with 



+ K +r/+T/)(Ap + — — ) 



(8) 



^^-ql XpFa - Va 



^ has exactly the same structure than r„ ^, but k is replaced hy Kf^ — n~ "E? jx- Note that in the absence of active 
noise, reduces to its thermal equilibrium expression, as it should. 

In the long wavelength limit, one finds for a membrane under tension (neglecting the osmotic contribution): 

{u{q^,t)u{-q_i_,t)) = — 2" (9) 

with T"^ = T (1 + IGApi^Q ?7^/xct) which means that a tense membrane is flat at longwavelength, even in the presence 
of non-thermal noise, and that one can define an effective temperature higher than the actual one. 

For a tensionless membrane, we now find for q± 0: 

{u{q^,t)u{-q^,t)) = ^ (10) 

nql 

with T''^ = T ■ {\pFa)'^K/[xiD + pXpF'^){p\pF'^ + ^XpFa)]- This expression is equivalent to the one given in Q|, and 
in particular the effective temperature does not depend on the pumping density for small enough densities. At long 
enough wavelength, the osmotic term should always dominate, but we show in the next section that the experimentally 
relevant regimes in fact imply the contribution of the force dipoles. 

IV. EXPERIMENTALLY RELEVANT REGIME AND DISCUSSION 

Let us first point out that in all equations the active-permeative and the active-hydrodynamic terms come in the 
combination: 

XpFa - -i-q^-w 

4:1] 

This tells us that in the long wavelength limit, the active-permeative term always dominates over the active- 
hydrodynamic term. However, the crossover wavevector, below which the active-permeative term wins reads: 

4.ivXp)Fa 

q±c — 

or for the corresponding length: 

That is with rjXp = Ip the permeation length, and Va/ Fa — w: 

2 

W 

l\r — StT 

Alp 

At first sight, one might be tempted to state that this crossover length is microscopic, but it turns out that Ip is of the 
order of, or smaller than, a Fermi. Indeed, with r/ — 10~^ kg/m.s and Ap ^ 10~^^m'^/N.s |Q , wc find Ip < 10^^^ 
m. Then, with w ~ 5 • 10^^ m, we find: 

~ 3 • 10^2 m 
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As a result, all active-permeation terms may be omitted in the micron and submicron length scales we are dealing 
with ^ 3 • 10^ m"! <q< 27r/a ~ 6 • 10^° m"!, with Ryes ^ 20 /zm and a ~ 0.1 nm). 

The osmotic contribution is negligible as well. To see this, compare 2XpkTa/ Dpq± and Paq±w/4:ri: this yields the 
cross over wavevector 

, _ rSXprikTa^ 

'^^^ ~ VaWDp 

or the cross over length 

VaWDp 



SXprjkTa 



1/2 



With an effective proton difFusivity Dp ~ 10^^ m^/s, a pumping rate d = 10"^ s^^, Va — n — lOfcT (see appendix |^), 
and other material parameters as above, one finds ~ 10~^ m: the osmotic contribution is totally negligible as well. 
In the case of ion channels for which d = 10^ s^"'^, the cross over length is reduced by a factor hundred, that is to a 
few tens of microns; this may be accessible to experiment. Similarly, since Ip q^^ , terms arising from permeative 
friction may entirely be omitted in the equations. We have now: 

i- {(Tqi_ + kq\) with k = k - pV^w 



4j) 



A further simplification can be obtained with the remark that in our experimental conditions D <^ {aK^^'^rj^^ (i.e., 

more precisely, D < 10^^(ctk)^^^77^^ with D ~ 10^^^ m^/s |2^). This means that one can further ignore the diffusion 
term in the {u{q±,t)u{—q±,t)) correlation function. Under such conditions, equation (|^) reduces to the expression: 

(u(q^, i)u(-qi, i)) ^ — ^ + ^ \ ^ (11) 

aqi + K^ql x(o- + 'f7l)('7 + k'"?!) 

where k"^ = k'^ — VaW^/x 

From this correlation function, we can calculate the relationship between the areal strain and the membrane tension: 

Aa = ao - a = — I — — ^ — ) In ( — ) (12) 

or 




Sttk (To 

with 



Aa = ao~a = — ln( — ) (13) 



fl+ VW-^Va-' \ 



The functional relation (|l^) is identical to the one holding in the equilibrium case (except that the 'temperature' is 
an effective nonequilibrium noise level) and is clearly compatible with the experimental results. The whole theory 
accounts well for the experimental observations if we note in addition the following four results: 

• the value of the bending modulus, measured in the passive case, is insensitive to the BR concentration and equal 
to that of pure phospholipidic membrane, 

• the effective temperature is about twice as large as the actual temperature, 

• the effective temperature is essentially independent of the BR surface concentration in the investigated domain. 
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• the reduced effective temperature difference {T'^f^ — T)/T decreases by a factor 3 when 25 % glycerol is added. 



The first observation is easily explained, as detailed in append ix |b| . In order to discuss the second and third observa- 
tions one must estimate the different terms entering equation (|14|). We provide details on these estimates in appendix 
p. We expect: 

\Va\ ^ K 

|S| ~ wkT 
X ^ kT/p ~ kTP 

where I is the average distance between BR molecules. We have been able to vary the concentration p over approxi- 
mately one order of magnitude, that is roughly I from w to 3w. With such estimates, and chosing the signs in such a 
way that the system is stable, we expect: 

T^ff ^ l + (^ + 2)^ 



T (1 + ^)(i+(2-M:)^) 
with 0.3 < w/l < 1 and knowing that k ~ lOfcT, we find: 

ye// 

1.7<^<2.3 

which is in very reasonable agreement with experiment. Of course, the numbers chosen above have some degree of 
arbitrariness, but one can change them appreciably while retaining a ratio T^ff/T of order 2. For instance, S may 
be set to zero, keeping other values unchanged, and one finds 2 < T^ff/T < 2.7, which is less satisfactory but not off 
scale. 

Let us now turn to the glycerol dependence. It has been measured that a 25 % glycerol addition to water reduces 
the pumping activity of the bacteriorhodopsin by a factor 2.5 |^^. It is thus clear that both Pa and V'^ (and perhaps 
S) have to be reduced by a factor 2.5; all other parameters are essentially unchanged, as shown by the experiments 
performed with red light. The same type of estimate as before give the expected reduction of {T'^ff — r)/T by a 
factor three. 



The net conclusion is thus that our analysis provides a satisfactory account of the experiment, although it is not 
able to pinpoint accurately values for S, and Va- A more accurate experiment should reveal that the effective 
temperature should depend on BR density in a non trivial way. At low density, the effective temperature should be 
essentially equal to the actual temperature; it should increase proportionally to the density at moderate densities; 
eventually, at larger densities, it could even decrease after going through a maximum. Together with an independent 
measurement of it should allow us to measure at least Va and V'a- Our current accuracy does not allow for such a 
detailed analysis. 

Thus, the proposed analysis gives a natural interpretation of the experimental data. It is one of those intriguing 
cases in which terms nominally subdominant in wavenumber provide by far the leading contribution. This pecularity 
is due to the very small value of the permeation coefficient: in the experimentally accessible domain, the membrane 
is practically impermeable and the effects due directly to the force exerted by the active centers on the membrane are 
neghgible. 



APPENDIX A: OSMOTIC PRESSURE DIFFERENCE 



Since BR selectively pumps protons, we just have to consider the osmotic pressure resulting from the three dimen- 
sional proton density n(r,t): 

n = fcTn(r,i) (Al) 

The protons dynamics is described as usual by conservation equations in each half space, above and below the 
membrane: 
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Fin 

^(r,t) + V-J„ = (A2) 
J„ = TiV - DpVn 

At the membrane, the coarsegrained proton flux J„ • z in the membrane normal direction z, is given by the BR active 
transport: 

J„ . z — 7/'(r, t) a 

where a is the pumping rate, and ip{r,t) is defined in the main text. For a macroscopically symmetric membrane, "0, 
n and V are 'smafl' fluctuatin g q uantities. So the convective term can be omitted, as a second order correction. Now 
to the hnear order, equation (|A2|) becomes in 'hybrid' Fourier space (with obvious notations): 



(A3) 



iujn{z, qj_, uj) + Dpq\n{z, q_^,uj) - Dp-^n{z, qj_, lj) = 

-Dp^^{z ^ 0+,qi,c^) = -Z?p|f (z = 0-,qj.,c^) = 7A(qx,c^)a 

The solution to this problem is straightforward. One finds: 

/ o+ \ t \ 2V'(q±,w)a 

n(z = 0+,qx,^)-n(z-0 , q±, ^) = ^ ^^^^/^ 

so that 

The typical frequency over which ip{q±,uj) varies is Dq^. Eventhough Dp is an effective diffusion coefficient renor- 
malized by the time the proton spends attached to the hydrazoic acid resulting from the conversion of s odiu m azide 
in solution one always has Dp ^ D, and thus the term uj/Dp may be safely omitted in equation (A4). Then, 
one can equivalently write: 

2fcrV;(q^,^)d 

OT(q_L,cj)~ (A5) 

Dpq±_ 



APPENDIX B: ORDERS OF MAGNITUDE OF THE THEORETICAL PARAMETERS 



The long wavelength effective membrane curvature modulus k!^ , is as shown in section [11 



"2 

= K-— (Bl) 

X 

It is easy to convince oneself that the coupling term 5 in fact depends on the pumping activity. Let us first consider 
the passive case, and call the corresponding coeffiecient Sp. A BR molecule with a given orientation may 'prefer' a 
given curvature sign for several different reasons. The first and most obvious one is linked to a putative wedge shape. 
In such a case, one expects: 

in which R is the 'radius' of the BR molecule and 9 the wedge angle. The experiments performed with red light tell 
us: 

^««: 
X 

That is: 

n2 kT ( I " 
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At the highest densities Z ~ -u; ~ i?, so with k = lOfcT the experiment requires 9 <C 1/3, which is obviously a 
'weak' requirement: the absence of up-down symmetry in BR requires the existence of a wedge, but inspection of the 
molecular structure suggests that it is very small (for instance, it is very hard to coin a sign to it). So, the 'steric' 
contribution to Sp can be safely neglected. 

There can however, be other contributions and the next most obvious one results from flexoelectricity: a curved 
membrane generates an electric polarisation, hence a transmembrane electric field. Again, the absence of up-down 
symmetry in BR tells us that it must have a non zero electric dipole. The energy of the dipole in this transmembrane 
field, provides the coupling between curvature and ip. With the usual definitions, the transmembrane electric field 
can be written: 

E = — -A^u (B2) 
ew 

in which e is the dielectric permittivity of the hydrophobic layer and e is the flexoelectric coefficient discussed by 
Petrov for instance [E^. If we call p the BR average longitudinal dipole, one has: 



5/ = 

The flexoelectric coefficient e is a measured quantity |2 



^ ew 



\e\ ~ 1.3-10"^" C 

Estimating p ~ a few q ■ S, in which g is a unit charge and S a distance of the order of a fraction of the membrane 
thickness, 6.17. the hydrophilic part (note that it cannot be much larger otherwise the BR would not be membrane 
soluble), and taking the dielectric permittivity of the hydrophobic layer of the order of e ~ 3eo, with the dielectric 
permittivity of vacuum eo — 8.8 • 10^^^ F.m^^, we find: 

=^ < 10-2 

So in this case as well, one does find k'^^-^ ~ k. 

When the BR undergoes its pumping activity, the flexoelectric energy is dominated by the time average energy of 
the proton in the flexoelectric potential eAj_u/e; assuming a duty ratio of one tenth, we expect then: 



96 1 ^2 



10 • e 100 • e 



q w kT 



The force dipole Va has the dimensions of an energy, and hence must be a fraction of the green-yellow photon 
energy. As a rough rule of thumb, we take again a duty ratio of a tenth: 

„ hv 
10 

The curvature dependence of the force dipole can be estimated in a way similar to that used for S. During its pumping 
cycle, the BR/proton system has probably to overcome a potential barrier Wb- The pumping rate is then controled by 
a Boltzmann factor exp(— VF^/fcT). In the presence of curvature, the barrier is modified by the energy of the proton 
in the flexoelectric potential at the barrier location. We call xWb this location. Thus, the activity is multiplied by a 
factor exp [—x{eqA±u)/{e ■ kT)] which can be linearized for small curvatures. Hence, we expect (with x of the order 
of a few tenths): 

Note that one could in principle estimate this coefficient by measuring the pumping activity in liposomes, as a function 
of liposome radius: the net result would strongly depend on the value of x. For x ~ a few tenths, one would need 
% accuracy to measure the curvature dependence. For x ~ 1, the effects would be much larger and the exponential 
nature of the relation should start to show up. However, even if the direct effect is not easily measurable, the incidence 
on formula dll) can be important. 
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